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The quantization of the Hamiltonian for a scalar field is performed in the framework of Quantum 
Reduced Loop Gravity. We outline how the regularization can be performed by using the analogous 
tools adopted in full Loop Quantum Gravity and the matrix elements of the resulting operator 
between basis states are analytic coefficients. These achievements open the way for a consistent 
analysis of the Quantum Gravity corrections to the classical dynamics of gravity in the presence of 
a scalar field in a cosmological setting. 


I. INTRODUCTION 

Quantum reduced loop gravity (QRLG) is a framework that describes the gravitational field of systems, whose 
spatial part of metric and whose dreibein’s are gauge-fixed to a diagonal form. It was introduced in mm and 
developed in m-m The theory has been successfully applied to an inhomogeneous extension of the Bianchi I model 

QRLG is constructed from Loop Quantum Gravity (LQG) [51I11| by imposing weakly gauge-fixing conditions in 
the kinematical Hilbert space. Therefore, it is a direct application of LQG that describes quantum cosmology and 
it differs from Loop Quantum Gosmology (LQG) |12II14] . in which quantization is performed in minisuperspace, i.e. 
after reducing the phase space on a classical level. The semiclassical limit of QRLG reproduces the quantum constraint 
adopted in LQG [3] and it also leads to enhanced inverse volume corrections dlH]. 

The most significant implication of LQG and QRLG is that the initial singularity is replaced by the Big Bounce m- 
The only possibility to test such a prediction is via the modifications it implies in the scalar spectrum of perturbations 
[I6l - [T8] . These modifications are due to the quantum corrections to the dynamics of the gravitational and infiaton 
fields. However, till now QRLG has been realized only in vacuum. Vacuum solutions in cosmology do not have real 
physical meaning for two reasons. First, Universe is filled with matter, which plays a peculiar role during its evolution. 
For instance, the photons of cosmic microwave background radiation provides the best source of information on its 
structure, while an inflationary phase can be realized through a slow-rolling scalar field. Second, time in General 
Relativity (GR) is not an observable and a natural way to account for it is via the introduction of a clock-matter field. 

In this paper, we introduce a scalar field in QRLG and we define the operator corresponding to its contribution 
to the scalar constraint. We quantize the field according with the LQG procedure given in [muni. Basic quantum 
variables are point holonomies and smeared momenta leading to polymer representation |211424] . The associated scalar 
constraint is quantized via a regularization of the classical expression, which provides a constraint written entirely 
in terms of SU{2) and point holonomies, together with the corresponding smeared momenta. This formulation is 
completely under control technically. This formulation is here adapted to QRLG, where the volume operator, thus all 
the relevant computations, are analytic. The final outcome of our analysis are precisely the analytic matrix elements 
of the scalar part of the scalar constraint between the basis elements of QRLG, which is the starting point for future 
applications. 

In particular, we give an introduction to QRLG in section |ll]by defining all the relevant structures of the kinematical 
Hilbert space. We focus our attention on states based at graphs having six-valent nodes |5] , which allow us to construct 
a cubulation of the whole spatial manifold. In section jlllj the classical and quantum formulation for a scalar field 
is given. On a classical level, we write the contributions of the scalar field to the scalar and vector constraints. 
The regularization of the field contribution to the scalar constraint is performed in section |IV| We just adapt the 
regularization performed in im to our case, which means replacing the triangulation with the cubulation of the spatial 
manifold and SU{2) group elements of LQG with the corresponding 47(1) group elements in QRLG. Having written 
the geometric variables in terms of fluxes and holonomies and the phase space coordinates of the scalar field in terms 
of point holonomies and smeared fluxes, the quantization is straightforward and it is performed in section |Vl The 
resulting operator is discussed in the large j'-hmit. We point out how if proper semiclassical states are constructed for 
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the scalar field, then the expectation value of the field contribution to the scalar constraint reproduces the classical 
expression. This result provides a first check on the consistency of the adopted framework. 

In this article we use the convention with metric signature (—, +, +, +), gravitational coupling constant k = IGttG 
and c = 1. Metric tensor is defined as g^i, = e^e^, where are vierbein fields. Dreibeins are denoted as e\, 
where lowercase latin indexes a, 6,.. = 1,2,3 label coordinate on each Cauchy hypersurface constructed by ADM 
decomposition |25], while = 1,2,3 are su(2) internal indexes. 


II. QUANTUM REDUCED LOOP GRAVITY 


The phase space of LQG is described by holonomies of Ashtekar-Barbero connections , smeared along some curve 
7, := V exp(/^ A^(7(s))r^ 7 “(s)^ and by fiuxes of densitized triads across some surface S, E{S) := UjtabcEjdx^A 

dx‘^. The kinematical Hilbert space of the theory is constructed as the direct sum of the space of cylindrical functions 
of connections along each graph T, 

»<liA=©«?’'’=i2(Adf..47. (1) 

r 

where A is the space of connections, d^AL denotes the Ashtekar-Lewandowski measure EH, while the states are 
cylindrical functions of all links S T and they are defined as 4'rj(A) := (A|r, /) := /(/lij(A), hi^{A), ..., for 

some continuous function / : SU{2)^ —C. 

The basis states labeled with a graph T, with irreducible representations D^’^ihi) (Wigner matrices) of spin j of the 
holonomy along each link I, and with an intertwiner implementing SU(2) invariance at each node v, are called spin 
network states and are given by the expression: 

= l[t, - (2) 

uGT i 


where the product U; extends over all the links I emanating from v and the • denotes contraction of the SU(2) indexes. 

QRLG implements the restriction to diagonal spatial metric tensor and triads along some fiducial directions, along 
which we define some coordinates x, y, z. The metric tensor reads 

dl^ = ttidx^ + a^dy^ + a^dz"^, (3) 


where the three scale factors are functions of time and of all spatial coordinates. The graph T now contains only three 
three kinds of links, each one being the set of links k along a fiducial direction. Inverse densitized triads are fixed to 
be diagonal. 


El=p^6l \p^\ = 


«iQ2a3 


( 4 ) 


(indexes are not summed in this expression), and this implies a SU(2) gauge-fixing condition in the internal space. 
Such a gauge-fixing is realized by the projection of SU(2) group elements, which are based at links k, onto U(l) group 
representations obtained by stabilizing the SU(2) group along the internal directions ui = Ui, with 


ui = (1,0,0) M2 = (0,1,0) U3 = (0,0,1). 


( 5 ) 


The kinematical Hilbert space now reads: 

( 6 ) 

r 

where T is a cuboidal graph, whiledenotes the reduced Hilbert space on a fixed (reduced) graph. The basis 
states are obtained by projecting SU(2) Wigner matrices on the state of maximum or minimum magnetic 

number mi = ijj, for the angular momentum component Ji = J ■ ui along the link h 

= {rniAi\D^'{hi)\rnuUi) , /i; G SU( 2 ). ( 7 ) 

Then, reduced states, called reduced spin network states, are given by the formula: 

= {h\{T,mi,iv}) = (jj,u;) • (/i;). rni = ±ji, 

uGT I 


( 8 ) 
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where (j;,uj) are reduced (one-dimensional) intertwiners. 

The graphical way to construct the elements out of those of the full theory is to replace SU(2) basis 

elements with the following objects 

\-(h^ 1 -^ = {ji,m\m'',ui) {m'',ui\D^’^{hi)\m'',ui) (m",m'), m" = ±ji. (9) 

Finally, reduction of canonical variables to and ^E{S) is obtained by smearing along links of reduced, cuboidal 
graph r and across surfaces S perpendicular to these links, respectively. 

The scalar constraint operator, neglecting the scalar curvature term, is obtained from that of LQG by considering 
only the euclidean part and replacing LQG operators with reduced ones. Its action on three-valent and six-valent 
nodes has been analyzed in [3] and [5], respectively. 

In what follows, we will consider the generic case in which the nodes of T are six-valent and we will represent the 
states based at each of them in the following graphical way: 


Jx,y,z 


\^-.U^)r = 



Ay) 

^ Jx,y,z 


( 10 ) 


where the node Vx,y,z is placed at (x, y, z), while the symbol j f. denotes the spin number attached to a link along 
the i-axis, beginning at the node {i + n,j,k} and ending at the node {i+n + 1, j, k}. We assumed the right-handed 
orientation of links, i.e. the link with the spin number j j. is outgoing from the node i n,j, k, while j j. 

is ingoing to the same node. 

In what follows, we will need the expression of the powers of the volume operator V, which acts diagonally |3] as 
follows: 


(y{vx,y,z)^ \T-U^)^ = {SttjI 


t .(x) 

Jx-1, 


y,z 


■ (x) 

■ Jx,y,z 


Ay) I Av) AA , AA ' 

Jx,y-l,z^ Jx,y,z Jx,y,z-1^ Jx,y,z 

X ■ ■ - -X 

2 2 


|r; u^l 


R > 


( 11 ) 


where 7 denotes the Immirzi parameter. 


III. LOOP FRAMEWORK FOR QUANTUM SCALAR FIELD 

The action of the scalar field minimally coupled to gravity reads: 

SW = 1 /■ d4xy^(5^"(5^0)(a.<(.) - Vi^)), (12) 

JM 

where A is the coupling constant of dimension 1/h and g is the determinant of four-dimensional metric tensor. 

The Legendre transform gives the following Hamiltonian: 
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N and N°‘ being the lapse function and the shift vector, respectively, while and Hsf'^ are the contributions of the 
scalar field to the vector and scalar constraints, q denotes the determinant of spatial metric and tt is the conjugate 
momentum to the scalar field. 

The total vector constraint is the sum of := Trda4> plus the gravitational part and it generates diffeomorphisms. 

The second term Hst\N] in the expression (131, which is the field contribution to the smeared scalar constraint 
encodes all information about the dynamics of the scalar field in the diffeomorphisms invariant phase-space. It can 
be written as follows 


:= / d^xN 


A 


Vd, 








r(0)r 


where we split it into three parts, the kinetic, derivative and potential ones, i.e. 

A 




(14) 

(15) 

(16) 
(17) 


We quantize the system of gravity and the scalar field by adapting the procedure described in [TOl EO] for LQG to 
the case of QRLG. Hence, the total Hilbert space is the direct product of that for gravity times that for (j) 


njdot) _ Rn/igr)^ ^{4>) 

’^kin ~ '^kin ^kin^ 

the latter being the following Hilbert space 

■“ + o,nU^^ : Oi G C, n € IN, G E,}, 

where the states are defined as 

with the normalization provided by the scalar product 


(18) 


(19) 


( 20 ) 


( 21 ) 


The basic variables act as follows: 


Ui, \U^') = \U ^+^,), n(H) \u^) = \U^) , 

vev 


( 22 ) 


n(H) being the scalar field momentum smeared over the volume H C E. One can also define single-point states 

\v;U^) := , (23) 


for which the scalar product reads 




(24) 


while basic operators act as follows: 




; t/y,) = |u; U^) = |u U w; U,p) , n(u) |u; U,p) = \ 1^’ 1^; > (25) 


d 


d(j){v) 


and one can define the smeared field around a point v via 

n(u) := jd^uxe{v,u)'n:{u). 


( 26 ) 







5 


where one introduces the characteristic function Xs{v,u) of the box Bg(v) centered in v with coordinate volume 
precisely 


V(B,(u)) := V(u,e) = e^y^iv) + 0(8^^), (27) 

which allows to smear a function at the point v, around infinitesimal neighborhood, such that 

f{v) = J(fiu6^{v-u)f{u) = \iin^Jd^ux6iv,u)f{u). (28) 

This way, the commutators are finite 

{(/)(a;),(/)( 2 /)} = {n(a;),n(y)} = 0, (29a) 

=X6{x,y). (29b) 


The full Hilbert space := L 2 (HBohr^) can be obtained from the single-point one T 2 (K,Bohr), where IR,Bohr 

denotes the Bohr compactification of a line and the Bohr measure is defined as 



rfMBohr(<?i’)e*'^' 


— ^0,1; • 


(30) 


This method for treating a scalar field on a lattice realizes a polymer representation in the momentum polarization 
(also called point-holonomy representation) |21fEl] . In QRLG, the only difference with respect to the scalar field 
quantization in full LQG is just the restriction to cuboidal lattices. 

The states in the total Hilbert space are described in the following way: 


|r ,Tfiijiyj — |r, TTii-f ^ |r, — 


,-(z) 


• ‘^'4^x,y,z+l4^x,y,z+l 



Av) 



A A 

Jx,y,z— 1 


o'^'^x,y,z-l4^x,y,z-l 


(31) 


The scalar field state is given by attaching at the each node Vp^g^r G T the point holonomy with the real 

( 7 ) (i) 

coefficient V'p.ij.r) while the gravity part is described by the spin numbers j] j at the associated links and the 
reduced interwiners at nodes. 


IV. REGULARIZATION OF SCALAR CONSTRAINT 


The quantization of the scalar part of the scalar constraint requires a regularization of the first two terms in the 
expression (141. This is done going back to the classical phase space and rewriting (141 in terms of holonomies and 


fluxes for the gravitational part, point holonomies and smeared momenta for the terms involving the scalar field. 

The gravitational part is regularized by the method developed in LQG [28], restricted to a cuboidal graph mm- 
Basically, the idea is just to replace the triangulation of the spatial metric with a cubulation. In this description, 
matter coupled to a dynamical spacetime is regularized by a reduction to matter fields coupled to a dynamical lattice 
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(nodes and links). Alternatively, in a dual picture, matter fields are coupled to dynamics of granulated space (volumes 
and areas of chunks of space) |8()| . 

In order to regularize and we mimic the procedure given in nanD] for full LQG in the presence of a 

scalar field. At first, smearing momenta and inserting the expression ^ = 1, where e is a determinant of triad one 


can write the following expression for the kinetic term (151 


Jd^xN{x)TT{x) Jd^yn{y) jdh 




{Y{u,e)Y 


-Xsix,y)xeit,x)xsiu,y), 


(32) 


where it has been used the definition of the volume V {R) of the region R 


V{R) = [d^x^. 
JR 


(33) 


Then, to remove the denominators in formula (321, we use Thiemann’s trick |28) . Using the Poisson brackets 


between connection and volume, one can derive the following relation 


eL(a;) = 2 


5Y{R) _ 2 5{Y{R)f _ 4 

n{Y{R)f~^ SEf nxK{Y{R)Y 


:i{Al(x),(V(i?)f} 


(34) 


and use it to write 


IdH- 


{Y{t,e)f 


— g| iijk 


1 

6\xk 


-A 


-A 


{Yit,e)Y {Y{t,e)Y (V(t,e))^ 


e^Jk J {A\t), (V(t, e)y } A { it), {Y{t, e)f } A { A'=(f), (V(t, e)f } , 


(35) 


which can be inserted into (32), thus giving 


A 2 


jd^r N{r)'x{r)xe{r,x) Jd^sTr{s)xe{y, s)x 

X J (V(t, e))^ }A{A^(t), (V(t, e))^ }A{A'=(t), (V(t, e))^ } x 

X J {a'(u), (V(u,£))^}a{a-(u), (V(u,e))^}A{A-(u), (V(u,e))^}x 
X Xe{r,s)xs{t,r)xe{u,s). 


(36) 


A method of discretization of the scalar constraint via a triangularization of the spatial manifold has been developed 
in |28| for pure gravity and then applied in the presence of a scalar field m- The idea is to replace the integration over 
the spatial hypersurface with the sum over over all ordered tetrahedra. Hence, the sum over tetrahedra becomes 
the sum over all the nodes v of the triangulation and over all the tetrahedra created by triples of links {Z, I”} 

emanating form v. Given a cubulation, each node v is always surrounded by three pairs of links, oriented along fixed 
perpendicular directions. They always create eight tetrahedra around the node and it is worth nothing that for each 
tetrahedron, the remaining seven ones coincide with the seven "virtual" tetrahedra, which must be constructed to 
triangulate any cuboidal or non-cuboidal lattice. 

Finally, integration over each tetrahedron, ^ , turns into the sum over the eight possibilities for choosing a 

triple of perpendicular links {I, I', I''} among each tetrahedron of the triangulation A(v) around the node v. 
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Using the triangulation procedure one gets the result 

(V(t,£:))"||A^(t), (V(t,e))" {V{t,e)Y^X6{t,w) = 

a{^,"(u), {Y{vi^i,,,,e)Y]sf„ <5„,„ « (37) 

~2"E E | (v(A(u), e))' , /ii^(A)|) | (v(A(z;), e))' , /ii,(A)|)x 

DGV(r)A(w) J/\L J/ 

where /i/(a) is the SU(2) holonomy operator in the fundamental representation |^, tr denotes the trace over SU(2) 
algebra and = —^a^, while are Pauli matrices. The summations X)iiGV(r) Sa(ii) extend over all nodes of 
the cubulation and over all the tetrahedra around each node, respectively. In the last line the following expansion 


tr(U/*-i{V"(i?),/r,J) =-tr(r*e{A,,V"(i7)} +0(^2)) « ^e{AlV-{R)} 
has been applied. As a result, the kinetic part of scalar Hamiltonian constraint operator reads 

221a 


tt( 4>) _ 
^kin 


E A(u)n(u)n(u')xe(i;,u') E E ^ijk^pqr^lmn^stu ^ 

['-fK) A'{v') 

^ {(v(A(u),e))', /ii,(A)|) |(v(A(i;),e))', /ii,(A)|)x 

X V(A)|^ tr^r'/i-J^,) |(v(U(A),e))', /ii,(A')|)x 

X I (v(U(A),e))', /izaa')}) { (vK(A),£))', /i/„(A') 


where (261 has been used. 

The derivative term (16) is regularized by the same method, applying the identity 

1 


/Xoafe _ _ 1 fc i Mf I m 


and it can be written 




{Vix,e)y {V{x,e)y 


X (fye)^e^^ddb(l3iy) 


{V{y,s)y (V(y,£)) 


-Xe{x,y), 


where the determinant of metric has been smeared and changed into a volume. 
Next, using again the Thiemann’s trick (341 one gets the formula: 


j pj 


€ijk(A I d-^xda(i){x) 


(V(j/,£)) 


r = f^ijk / d 4 >{x) A 


-A 


16 

SyK 


{Y{x,e)f {Y[x,e)Y 
ez3k fd(l){x)A^A\x), (V(a:,e))"|A|A'"(x), (V(a;,e))"|, 


which applied to the expression (411 leads to the following result: 

rrW ^ 1 216 


dev - 8A 34(7K)V bfc6’;m 1™ jN{x) d(j){x)A^A^ [x), (V(x, e))^ } A |a'=(x), (V(a:,e))^|x 
X fd(j){y)A\^A^y), (V(2/,£))" |A|A'=(y), {V{y,e)Y'^Xeix,y). 


(38) 


(39) 


(40) 


(41) 


(42) 


( 43 ) 















By using the same discretization adopted for the kinetic term (371, one gets 


(V(x, e))'‘ {V{x,e))*'^Xe{x,y) = 

E ^ijk^pqr <9j3 0('u) tr 

uGV(r) A(t;) 



(44) 

j 


and for ff 


^ 34A hnY i™ E E E ^ijkipqrAm^^stu 5v,v' X 

^ ^ ’ «,'u'GV(r) A(j;) A'(v') 

X 5p(/)(u)tr(^r^h-E |(v(A(n),e))", /i/,(A)|^ tr | (v(A(n), e))" , V(a)|^x 

X ds(j){v')tx{T^h-\^^ |(v(A'(^;'),e))^ | (v(A'(r;'), e))^ /iz„(A)|^- 


(45) 


If the scalar field is sufficiently smooth one can write 

1 e'^”+=p"‘^’' - 


dp(i){y) 


(46) 


where (pv+e is the field in the point v + e),, which is the nearest node of v along the link of length e. Finally, 


applying it to the formula (45) one gets 

ol3 

tW _ ^ 


- 34A(^k)4 E E <^iikf-pqri\m^stuy- 


i’GV(r) A(4;)=A'(t;)=i; 
ct 


v + ^p-<Pv _ ^(pv-4>v-ep + 


X 


X {(V(A(^;),£))^ /i/,(A)}) tr(^''V(A) {(V(A(u),£))^ / 14 (A)}) x 

X {(V(A(^;),£))^ / 14 (A)}) tr(^'”/i);(A) {(V(A(u),£))*, /i 4 (A')})- 


(47) 


V. QUANTIZATION OF THE SCALAR HAMILTONIAN CONSTRAINT 


The field contribution to the scalar constraint is quantized by the canonical procedure: the cubulation of the spatial 
manifold is given by the graph F at which the state is based (links and nodes of the cubulations are links and nodes 
of F), while holonomies, volumes and matter variables are changed into quantum operators that act on states (311 
belonging to 


H\T] rriiAv] “ {^ktn + + -^pot) 1^’ . 


( 0 ) 


r( 0 )) 


(48) 


The Poisson brackets in (39) and (47) are replaced by commutators and the following condition is going to be used 




v"(i?), K 


= tr[T^hiy^{R) \ 


(49) 
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The quantum operator corresponding to the kinetic part (391 acts as follows: 

221a 




^ 32(167r7C?/i)® e->-o ■■ 

1 ' ' «.j)'GV(r) 


lim E Ai(u)n(u)n(u')x£(i'.w') E E ^ijk^pqr^lmn^stu ^ 

A(if) A'(v') 


X til (V(A(u), £))' /r,,(A)j X 

tn ('^'"KIa) (A(t’), ff))" hiAA^ tr (v(w'(A'), e))" K{A'^ x 


(50) 


X tr (v(u'(A)> £))' ^/t(A')) tr (i'"^„(A') KA). e))" ^i„(A')] jT; mp t/^)^ , 


where V is the volume operator (111 and n(iy) is the momentum operator, whose action in Schrodinger representation 
reads 


n(i(;) \r;mi,iy; = -ih 


d 


d4>(w) 


|r, TTLl , 2^, . 


(51) 


H, 


( 0 ) 

kin 


2^0^ _, „ _, 

|r, 771/, S^TStT \6 ^ / A^(u)n (u) ^ ) ^ijk^pqr^lmn^stu ' 

t ' Pl ..cvrr'i A7„'i=A'0,'i=,) 


Then, taking the limit e —>■ 0, in the expression (501, one reach the scale at which there remains only the single 
node V = v' and the dependency on the regulator e is removed: 

2i5A 

,V 

j;GV(r) A(v)=A'(v)=v 

X tll(^*E(A)('^(^))'E(A)) ^^(:^'K}A)i'^(^)yK(A)) X (52) 

X ^^(;^'"K\A)(y(^)YKiA^ tl'(l‘*^/”(A')(V(p))"^iAA')) X 
X (V(p))" ^q(A')) tr(r”/7/;'(A) {V{v))~^ KiA)) |r; . 

The result is a sum of subsystems, called basic cells, that extend over all nodes of the graph T. Each basic cell is 
a sum of elements acting on each node surrounded by six nearest neighbor nodes and it is labeled by the position of 
the central one. The example, with the central node Vx^y^z-, is given by the illustration of state (31). This cellular 
structure allows to restrict calculations to a basic cell and to give the final result as the sum over ml cells. 

The action of the operator (52) can be computed from the following expression (see appendix) 


tr(r*/7-iV"(u)/i/J|r;777/,7„;[/,^)^ ='- S^p , 


(53) 


where we introduced 




^ [(I# - 1 |++ 1 | + 

. . . . .(p) 

Cp being the unit vector along the direction p, such that j]^_g is the spin number of the link along p ending in v. 


(54) 


Hence, the kinetic operator (521 reads: 


777/, 2„; C/^)^ — 


23 a 


32 ( 8777 / 2 ,) S 

E E ^ 

A(v) {p,< 7 ,r}GA(4;) 


^fv.ri2 s^) s(^)x 


{p),\ A (9),i A W.i 

L\v AAd 




(55) 


A'(i;) {s,t,n}GA'{t/) 


The summation J2{pqr} comes from j k p q r^ijk^pqr^ip^jqSkr in ( [5^ (with Kronecker S’s arising from (53)) and 
extends over the 6 permutations of links in a given triple. The summation X)a(j;) ^cie to the 8 possible choices 
of triples with mutual orthogonal links, which span the 8 tetrahedra A surrounding the six-valent node. Since 
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(p) - (<?) - ('T') - 

Ay A„ ^ is invariant under permutations of links and choice of tetrahedra, the final expression is just 

(6 X 8) that of a single triple. In the same way, we treat X]a'(d) X]{s t u}gA'(?;)- Finally, we get (6 x 8)^ times the 
eigenvalue of the product of three trace operators acting on a single triple: 


mi,Zy- 3 Y^Ny fig 

( StT^/^ I 2 - \ 


{iTT-iPp)2 


i A(y).i 

Z-it; i-Av 


|r, uii , 2^, ■ 


Similarly, the quantization of the derivative part of the scalar constraint operator (471 gives: 

213 

t;eV{r) A(i?)—Aq-i;)—f 

d 

— 


-^der|F; nr/, li,; ^4 ^ijk^pqr^lra^stu^^ 


v-\-e.p 


X 


ep g9. + e 

— X - 


— 


-X 




Next, using the expression for the trace ( |53[ ) one obtains: 

^ ]V(u)^ 51 51 51 e^Jkepqre^m^stu6qJSrkSltSy 

tiGV(r) A'(f) {s,i,ii}GA^(t;) 


t,+Sp-4>^ _ p4>-v-4‘-a-s, 


■p f>9p+P 


(56) 


(57) 


(58) 


(sw)^ 


aI’'^’® aI*^’® Ai“^’®ir: 




|r, TTll , iy , U-lp)^ . 


Since the expression within the summations is invariant under the exchange of q i—> r (and of t i—> u), we get 

H^/2\T; mM.; U^)n = J2^v {^2^ ^2^ J2 ^ 


A' (v),A{v) 


^0u + e„-<^' 


p 'T’v _ ^yJv—iPv-ep ^(pv + es—<Pv _ 

-X --- 


(59) 


(EWE^"i)'(Ai'^^’®Al’'^’® 


3\2 


|r, 77T./, , 


where 6ps arises because the J’s in (581 force p = i and s = i. Note that the number of terms involved in the 
summations differs from the one in the kinetic part of the Hamiltonian. The summation extends over eight tetrahedra 
in the both cases, A and A, giving 8^ terms (as in (55)). However, the Kronecker delta <5^^ identifies one edge of A 
(Ip) with one edge of A' (4). Therefore, the action of Hamiltonian extends over eight tetrahedra A'(u) and then over 
the four tetrahedra A(u), which share the link Ip. 

Finally, contracted indexes p and s give the summation of over three directions {x,y,z} and the equation (59) 
becomes: 


Sllr; [/,)„ = (£(,-)£(»)£(•))• x 

-T("F)*( 


+ 




ey — g^u + ey 




2 

ex _ + ‘ 


(£<»>)'( 

( 4 '')'( 


^y ^y 


(60) 


A(^)’ 8 A^^^’ 8 ^ 


|r, 777./, iy , P'tp) 


R ■ 
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The same result can be also obtained in other way, directly form the formula (571, by the following counting: there 
are six possibilities for attaching a first triple of spins to the links {lp,lq,lr} C A(n) times two possibilities for the 
second triple of links {ls,hjlu} C A'(n), while the summation over triangularizations extends over eight tetrahedra 
A'(n) and over four tetrahedra A(n) (which have the link 4 as one of their edges). 

As for the potential part quantization is straightforward, since the only geometric part is the volume V(n). 
Using (111 one gets the expression: 




«GV(r) 


(SttjI 


2 


2A 


|r, nz/, j 


(61) 


where the operator V(cj)) should be properly defined in the polymer representation. 
Then the action of the total scalar constraint operator (481 reads: 




2“A 


2 'll 


yM yiv) v( 2 ) ( 

V V \ 


A A). 3 A(y)U a(z) 

AA 7 ; AAii AA 7 ; 


1 \2 


ni 


(4„)5(y,i,l)»(^(4l)» 






-^y — q^v + E 


(62) 






(STryZ 


2 'll 

p) 


(si") si^) sl^))" IT; mi, ly -. 


The expression above gives the action of on a quantum level and it is the starting point for the analysis of the 
dynamics of the scalar field. It is worth noting that all the coefficients within (621 are analytic. 


V.l. Large j limit 


Let us now perform the large j limit of the formula (62) and outline how the eigenvalue of the quantum Hamiltonian 
coincides with the classical expression (14) at the leading order. To calculate this limit one can consider the following 
expansion for j ^ 


AW’" = -n(sW)" Vo(j—3), 

and gets for the expectation value := {T;mi,iy-, U^jj] H^^^\r-,mi,iy; 


(63) 




Snjll Sl^) y 


2 vU) vU) 




SiTjll si") si^) 




^Vy {V{M)\ 


(64) 


where we introduced the eigenvalues V^, and n„ of the volume operator V and of the momentum operator H^ 


V„ := 


:= {^{STT^llf e( 9) El"))" n„ = hiPy 


(65) 
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and the expectation values (Af(^) and i.e. 


= {U^\ 


«-Sp _ pVv + g 


m) {vm = {U^\V{cl,)\U^) . 


( 66 ) 


Next, using the definition of the characteristic function p8|), one gets 
1 


K, lim / d^u Xe(v, u) N{v) ( — ^ ^ 


+ e 




p^{v) 




p^(v) 




p^{v) 




+ e- 




(67) 


2A 




2A [p2(^;)p3(^,)' p^v)p^iv)^ p^[v)p'^{v) 

where {u) denote gravitational momenta at the point u, which are related to spin-numbers by the following relation 

p\v) , (68) 

and q = |p^p^p^| is metric determinant. 

Let us assume to construct a proper semiclassical state for the scalar field variabletQ such that expectation values 
and eigenvalues become classical quantities. Hence, the semiclassical value of becomes, in terms of the 
original un-smeared variables. 


iT-if fd^^Xe{v,u)N{v)(—^=TT^{v) + 


\/gM 

2A 


P^iv) /„ ,, ^^2 P^iv) /„ ,, ,s2 P^M t r. ,, ^\2 

n, , p, J dx(^{v)) + 1/ N A dyH'v)) + 

p^{v)p^[v)^ ' p^{v)p^[v)^ ' p^{p})p^(v)^ ' 




(69) 


and in the limit e —> 0 we have v = u and ^^Jd^uXei'^j ~ f finding 

Jd^uN{u) ^7r2(u)-H|^^9iya„(/)(u)f-hg2yaj;((>(u)f, 

where we introduced the inverse components of the metric tensor in terms of p’s 




p2p3 




p3pl 


33 P 

9 = — ■ 

pLpZ 


(70) 


(71) 


The expression above clearly coincides with the classical expression (141 with the metric in the diagonal gauge. 


VI. CONCLUSIONS 

We defined the action of the scalar field Hamiltonian in the diffeomorphisms invariant Hilbert space of QRLG, whose 
quantum states are based at cuboidal graphs with attached t7(l) group elements. Hence, we adapt the procedure 
defined in m to a cubulation of the spatial metric and to the reduced holonomies and fluxes proper of QRLG. The 
scalar field was described in terms of point-holonomies, which live only at nodes of the graph. The resulting action 
of the scalar field Hamiltonian has been regularized via standard tools of the full theory and its matrix elements 
are analytic. This is a key point of the formulation, which outlines how the dynamic analysis can be carried out 
analytically in QRLG not only in vacuum [3], but also in the presence of a scalar field. We expect that this is the 
case also for other matter fields, whose introduction in QRLG will be the subject of future developments. 

We also checked how in the large j limit, the Hamiltonian eigenvalues approach the classical Hamiltonian at the 
leading order. The next-to-the-leading order corrections are pure quantum corrections, which will be discussed in 
conjunction with their possible phenomenological implications. 


^ The construction of semiclassical states for a quantum geometry in the presence of a scalar field is in preparation, along the lines of 
what has been done in vacuum [3- 
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In particular, there are two kinds of applications of QRLG in the presence of a scalar field. The first one is to 
regard the scalar field as a clock-like field, defining the evolution of geometric degrees of freedom on a quantum level. 
This is analogous to what has been done in LQC min]. However, the quantization procedure described here is not 
equivalent to the canonical one adopted in LQC, thus we expect nontrivial results from this analysis. 

The second application concerns regarding the scalar field as an actual matter component of the thermal bath 
and to analyze the corrections to the classical dynamics. This investigation could be relevant for those inflationary 
scenarios, which are based on a scalar infiaton, since it can predict modifications to both the background expansion 
of the Universe and the behavior of scalar perturbations. 


Appendix A 


In this appendix we compute the action of the operator tr 
simplicity p = z, we have 




(v) hi ^) and we infer (531. Let us choose for 


= -Y,{n)ab{hT")bdV^ {hljda (Al) 

abd 

with a,b,d indexes in the fundamental representation and r SU(2) basis elements; in the basis that diagonalizes r^, 
the holonomies are diagonal: 

(A2) 

When we apply this object to a state, the volume acts after the insertion of the holonomy hi^, so it gives a coefficient 
sC) (s(U + a)]"/2 and we get 


in)abCK")bd +ay'UhiJaa 

abd—— 1/2 
1/2 


1/2 


n/2 




V2 




abd—— 1/2 

Now using the 6’s we get that a = d = b, such that the two exponentials disappear and 

n/2 _^ f t \ \ 


5](a:) j](y) 




(A3) 


(A4) 


which, by considering that the only r with nonvanishing diagonal components is r^, becomes 


% 

2 1 


y}'V) 


nl2 




n/2 


4 1 


y}'V) 


n/2 


+ 1/2)”^^ - -1/2) 


n/2 


■ (A5) 


From the expression above equation (531 for p = z follows. For p = x,y, hi is diagonal modulo the rotations, which 
can be moved to t^. As a consequence, we get the same result but with the rotated r, which means that the only 
non vanishing contribution is for i = x,y. 
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